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We introduce a simple model for DNA-cationic-lipid complexes in which galleries between planar 
bilayer lipid lamellae contain DNA 2D smectic lattices that couple orientationally and positionally 
to lattices in neighboring galleries. We identify a new equilibrium phase in which there are long- 
range orientational but not positional correlations between DNA lattices. We discuss properties of 
this new phase such as its x-ray structure factor S{r), which exhibits unusual exp(— const. In'^ |rl) 
behavior as a function of in-plane separation r. 
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DNA is a remarkable polymer that exhibits a complex 
phase behavior as a function of packing density, salt con- 
centration, and other variables [y . It is anionic, giving up 
positive counter ions to solution. Mixtures of DNA and 
cationic and neutral lipids in water form complexes that 
facilitate transfection of DNA into living cells and that 
play an important role in the emerging field of gene ther- 
apy H . Recent x-ray experiments reveal the structure 
of these complexes at length scales from 10 to several hun- 
dred Angstroms, particularly near the isoelectric point 
where the total charge of counter ions given up by the 
DNA equals that given up by the cationic lipids. The 
lipids form bilayer membranes that stack in a lamellar 
structure (Fig. |l|). Parallel strands of DNA arrange in 2D 
smectic structures in the galleries between lipid bilayers. 
The distance between lipid bilayers is equal to the diam- 
eter of a DNA molecule plus a hydration layer. In addi- 
tion, the distance d between DNA strands increases with 
increasing concentration of neutral lipids in a manner 
consistent with counter ions being expelled to solution 
and charge neutrality of the complex being determined 
only by the DNA and cationic lipids. The best fit to 
x-ray diffraction data is obtained when some correlation 
between DNA lattices in different galleries is introduced. 

We undertake here a theoretical investigation of possi- 
ble equilibrium phases of these lamellar DNA-lipid com- 
plexes. We identify a new phase, with a nonvanishing 
smectic compression modulus B in each gallery, in which 
there is long-range orientational but not positional cor- 
relation between DNA lattices in different galleries. This 
phase exhibits no restoring force for sliding DNA lattices 
rigidly relative to each other, but it does exhibit a restor- 
ing force preventing their relative rotation. We will refer 
to it as a sliding columnar phase. It is distinct from both 
the columnar phase in which the DNA segments form 
a 2D lattice and the totally decoupled phase in which 
there is no communication between different DNA lat- 
tices. It is similar to the decoupled phase of stacks of 
tethered membranes [0. Dislocations can destroy po- 
sitional correlations within the DNA lattices, melt the 
sliding columnar phase, and produce a nematic lamellar 
phase with B = 0. Whether they always melt the phase 
will be discussed in detail elsewhere H . 




FIG. 1. Schematic representation of DNA-cationic-lipid 
complex. Parallel strands of DNA form smectic lattices with 
lattice spacing d in galleries between lipid bilayers with spac- 
ing a. Charged and neutral lipid heads are, respectively, 
shaded and unshaded. DNA strands are aligned parallel to 
the X axis, and the y axis is the normal to the lipid planes. 

We consider a model in which the DNA strands are 
confined to galleries between lipid bilayers in a perfect 
lamellar structure (with layer spacing a) with no dislo- 
cations or other defects. We assume the ground state of 
DNA strands in each gallery n is that favored by elec- 
trostatic interactions, i.e., a 2D smectic lattice with layer 
spacing d — 27r/fco- We take the lipid bilayers to be paral- 
lel to the xz plane and the DNA strands to be aligned, on 
average, parallel to the x axis as shown in Fig. |l|. For the 
moment, we assume that the lipid bilayers are perfectly 
flat and do not fluctuate. In this case, long-wavelength 
properties of the DNA lattice in gallery n are described 
entirely in terms of displacements u" (r) along the z di- 
rection, where r = (x, z) is a position in the xz plane. 
The Landau-Ginzburg- Wilson Hamiltonian for the com- 
plex is then a sum of independent elastic energies for each 
gallery and terms coupling displacements and angles in 
neighboring galleries: Ti. = H''^ + E„(^« + ^n) with 

ni = -Ve f d^r cos[2(r - 0"+i)] (1) 



ni = -Vu / Srcos[kfi{i 



,n+l\ 



where 61" w 9^u", and <^ = a^'u"-[(92;M")^ + (92M")^]/2 
is the nonUnear strain for gallery n. B2 and _ftr2 are, re- 
spectively, the 2D compression and bending moduli. Vu 
is of order {\^/d)e~'^^°'^'^ where Ac is the charge per unit 
length (e per 1.7 A) of DNA. We do not yet have an es- 
timate of Vg whose dominant origin is likely a membrane 
mediated interaction. 

If Vu is sufficiently strong, the DNA strands form 
a regular 2D lattice in the yz plane, and the whole 
complex will resemble an Abrikosov flux lattice in a 
high-Tc material with the magnetic field parallel to the 
copper-oxide planes [|6| . It will have the symmetry of an 
anisotropic discotic columnar liquid crystal 0. As T4 is 
reduced (or temperature is increased), fiuctuations can 
melt the columnar DNA lattice without destroying the 
lipid-bilayer lamellar lattice. What is the nature of the 
melted phase? To answer this question, we consider the 
limit in which the potentials Vu and Ve are small. From 
the known statistical properties of 2D smectics, we find 
that Vu is irrelevant and Vg is relevant with respect to 
the totally decoupled phase. 

When the potentials Vg and Vu are zero, we have a 
stack of decoupled 2D smectics, whose properties are 
by now well known [plpl. First we will ignore disloca- 
tions. At length scales less than the nonlinear lengths 
l^ = Kl^'^/Ty/B^ and h = ll/\ where T is the tem- 
perature and A — \/ K2/ 82, fluctuations are described 
by the linearized elastic Hamiltonian with the nonlin- 
ear strain replaced by the linear strain dzU. At lengths 
scales longer than l^ and Z^, nonlinearities lead to renor- 
malized bending and compression moduli K{c[} and -B(q) 
that, respectively, diverge and vanish at small wavenum- 
ber q. In both the harmonic and nonlinear regimes, the 
displacement correlation function in each gallery can be 
expressed as 



G(q) - q-^Q{qJqi) 



q^' 



'^'^^ (2^ 



In the harmonic regime, 77 = 4, /i = 2, and TG^^{q) = 
^2(0^ + ^^Q-t)- I'^ the anharmonic regime, ry = 7/2, and 
/i = 3/2. The mean-square fluctuation in the displace- 
ment diverges in both regimes with the lengths L^ and 
Lz of the sample in the xz-plane: 



_fq_ 

(2^)2 



G(q)=Lr/W(i./LS), 



(3) 



where 2a = r/ — 1 — /i = 1 in both regimes and /„ (0) ~ 
T/^/K^ and f'i,^\w) - A^/J^^^^w^"/'^ as w; -> 00. 
Alternatively the displacement correlation function, 

gu{v) = {[{u{v)^u{Q)Y) - I^P"/?Hl^l/kr), (4) 

diverges with separation r. Angular fluctuations are 
non divergent : {9^) = / ^^^^(q) = A2"/,(AJA^), 



where A^, ^ \jd ^^ A^ ~ l/d are the wavenumber cutoffs 
in g^; and q^ and /^(O) = const, and Jq(vS) ^ it;^"/^ . 

Dislocations deflne another length scale ^20 = n^ ^ 
exp{E2iy/2T), where n^ is the density of dislocations and 
^2D is ^^6 energy of a dislocation, which is finite for 
smectics. At length scales less than ^^D' ^^^ system is 
described by the harmonic or the nonlinear elastic the- 
ory. At length scales longer than ^jdj dislocations melt 
the smectic lattice, leaving a 2D nematic with power-law 
angular correlations [0. In what follows, we concentrate 
on the case in which ^2d is the longest length scale in the 
problem and in which there is a regime in which disloca- 
tions can be ignored. 

The coupling energies 7i^ and TiJJ are irrelevant if they 
tend to zero at large L^ and L^ and relevant if they di- 
verge with Lx and L^- If DNA lattices are totally de- 
coupled, (H") = —VuLxLzexp{—k^{u'^{r))) tends expo- 
nentially to zero for both harmonic and anharmonic 2D 
elasticities, and is irrelevant. On the other hand (6'^(r)) is 
finite for both linear and nonlinear elasticity, and the po- 
tential (Tin) = -~VgLxLze~'^^^ ^ diverges with L^Lz and 
is relevant. When dislocations are allowed, then (cos0„) 
dies off as a power law at large L^ and Lz, and (7i^) 
may be relevant or irrelevant. Thus, if we ignore disloca- 
tions, the angular coupling is relevant, and the decoupled 
Hamiltonian will flow to a new long-wavelength Hamil- 
tonian with angular but not positional coupling between 
layers. In the continuum limit, this Hamiltonian is 



H = i f d^x[Bulz + K{dluf + Ky{dydxu)\ 



(5) 



where x = {v,y), u(x) = ^^'/"(r), B = Ba/a, K ^ iCa/a, 
and the nonlinear strain is still that of a two-dimensional 
system. Note that there are no terms proportional to 
{dyuf' or [dyuf' . The only term involving y derivatives is 
the one Ky{dydxu)^ = Ky{dy9)^ coupling angles in differ- 
ent layers. This Hamiltonian is invariant under transfor- 
mations of the type u(x) —> u(x) -I- f{y) for any function 
f{y) independent of r, i.e., it is invariant under arbitrary 
rigid displacements of one layer relative to another. To 
lowest order in Ve, Ky = Vga{cos26)l — Vgae^'^'^'^'K 

The elastic constant Ky introduces a new length ly = 
ayjK/ Ky. At length scales within a gallery less than ly, 
DNA lattices behave like independent 2D smectics. If 
ly < 'x,z) there will be a crossover from 2D harmonic to 
37? sliding lattice at length scales of order ly. If ly > l^^z, 
there will be a crossover first to 2D nonlinear behavior 
and then to 3D sliding behavior. 

Fluctuations in u in the harmonic limit are now de- 
termined by G'(q) = T[Bql + Kq^ + Kyqlql]~^. They 
diverge with system size L: 



C\n L and ((Awp 



ApC\nL, 



(6) 



where Ai 



— ,,n+P 



(r) 



'(r), C = Tl^JBKy, and 



Ap is a number. Thus, {JV^ w -Ki^"''% where 



rjc — AifcgC/2 tends to zero at large L for rjc > 2, and 
K is irrelevant for T > Tu = {'i/AikD^BKy. 

Since {v?) diverges with L, the correlation function 
5„(x) = ([u(x) — ?i(0)]^) diverges with x. In the limit 
L -^ cx), (7u(r, 0) is finite for all r and at large |r| is 



5„(r,0)~Cln2(A2A|r|). 



(7) 



When y is nonzero, 5„(r, y) diverges with L for all y. 
The behavior of (7«(x) has interesting consequences for 
the x-ray structure factor iS'(q) ~ J d'^xS{x.)e^'^'^'^, where 
5(x) = (eifco[«(x)-«(o)]<)^ Pqj. y ^ 0, the dominant contri- 
bution to 5(r) = S{r,0) = exp[-fcgg„(r,0)/2] is 



S{r) 



,-Cfe^(ln-' A^A|r|)/2 



|r| > I 
e-fcobl/f'(|2|/bl")/2^ |r| < ly 



v 



(8) 



The contribution to S{0,y) arising from guiO^y) is zero 
in the infinite volume limit when Vu is irrelevant. There 
are, however, short-range positional correlations between 
layers arising from the irrelevant variable Vu- To low- 
est order in a perturbation expansion in V^/T, 5(0, y = 
no) ^ (T/u/2T)"/" where / has contributions of the form 
/ (i^re^'^"^"'-'''^^/^. Thus, there is exponential decay of 
5(0, y) - e-\y\'^y with ^^ = a/ln(2T/K/)- The cor- 
relation length £^y is finite so long as Vul < 2T. Since 
.g„(r, 0) grows more rapidly with r in the totally decou- 
pled phase than in the sliding phase, / and ^y are smaller 
in the latter than in the former phase. There are further 
contributions to 5(r, y) that become more important as / 
grows larger. When ^j^ diverges, there will be a transition 
to the true columnar crystal phase with a nonvanishing 
2D shear modulus. 

The Hamiltonian of Eq. (|5|) has nonlinear parts arising 
from the nonlinear strain Uzz ■ These lead to logarithmic 
renormalizations of the coefficients B, K, and Ky similar 
to those in smectic liquid crystals jl^]. We find 

if,(q) ^ K^/^q) ^ i3-l/3(q) ^ [ln(A//i(q))]l/4, (9) 

where /i(q) = [q^ + X'^q^ + Xlqlql]^/^ with A^ = K/B and 
Xy = Ky/B. Preservation of rotational invariance in the 
calculation of the above renormalized elastic constants 
requires considerable care. Details of this calculation will 
be presented in a separate publication |l^]. 

We have thus far ignored dislocations. In two- 
dimensions, dislocations cause the smectic lattice to melt 
at length scales larger that ^^ . Dislocations will certainly 
suppress order in the sliding phase. To obtain quantitave 
estimates of their effect, we need to calculate the energy 
per unit length of a dislocation predicted by the elas- 
tic energy of the sliding phase. We would like to obtain 
not only the elastic energy of a dislocation but also some 
estimate of its core energy, which can arise either from 
destruction of the order parameter "0, describing peri- 
odic order, over an order-parameter coherence length f 
or from deviations of the local director (specifying in this 



case the normal to the DNA strands) from the direction 
of preferred alignment over a bend or twist penetration 
depth Xb- In the smectic- A phase, screw dislocations 
have no clastic energy; their energy is totally in the core 
and is dominated in type II systems, for which Xb > £,, 
by director misalignment |12|] . The mass-density ampli- 
tude ?Adna for our DNA lattices is strictly speaking un- 
defined in the sliding phase in the regions between DNA 
galleries. Thus, we argue that any core energy arising 
from dislocation lines parallel to the layers (i.e., in the 
xz plane) must come from director mismatch rather than 
from destruction of periodic order. To describe director 
mismatch, we consider a "gauge" version of the elastic 
energy of the sliding phase in which displacements and 
angles are coupled in a rotationally invariant way: 



Tin 



(fx[B{dzuf + D{d.^u- ef 



-K^idjf + KyidyOf + KAdzOy 



(10) 



This energy reduces to the sliding phase elastic energy 
[Eq. (w] when 9 is integrated out. It introduces twist 
and bend penetrations depths Aj, — \J KyjD and Xz = 
^J Kzl D characterizing length scales over which orienta- 
tions relax in response to layer distortions. 

Dislocations are topological hne defects. They can 
be characterized by a density b(x) that determines the 
singular part of the displacement variable u(x): V x 
Vu(x) = b(x). Using standard techniques ||l^, we 
can calculate dislocation energies in terms of the Fourier 
transform b(q) of b(x): 



^^2^ 



d^c 



K{q)q^\by{ci)\' 



{27:fBql+K{ci)q^[ql + {BlD)qlY' 



(11) 



where K{c\}q^ — K^q^ + Kyq^ + Kzql- This energy de- 
pends only on the y component of b, i.e., on the com- 
ponent of b normal to the bilayer planes. When there 
is a single dislocation line parallel to the y axis with 
b(x) = deyS{x)S{z), the energy per unit length e pre- 
dicted by this equation is precisely E2BI0. where E'ju is 
the energy of a two-dimensional edge dislocation. Since 
Eq. (O) includes orientational mismatch energy and we 
expect no core energy from destruction of the mass- 
density order parameter, we conclude that the energy 
for dislocations with b in the xz plane is in fact zero. 
Thus, it is possible to form dislocation loops with arbi- 
trarily long segments in the xz plane. (This situation is 
analogous to that of the "pancake" models of flux lat- 
tices in high- Tc superconductors [H). These loops are 
equivalent to independent 2D dislocations connected by 
loop segments with by — passing between DNA gal- 
leries. Dislocations will melt the sliding columnar phase 
if these 2D dislocations melt the smectic lattice in any 
given gallery. Smectic lattices in galleries above and be- 
low a given gallery act as external fields orienting that 
gallery along the common direction of the whole sample. 
Thus, each gallery is equivalent to a smectic in an exter- 
nal aligning field (whose long-wavelength Hamiltonian is 



the in the universahty class of the xy model [Q ) that will 
Kosterlitz-Thouless melt to a nem.atic at a tem.perature 
Tjv. If Tn > Tu, then for T^ <T <Tn, there lattices 
have not melted, Vu is irrelevant, and the sliding colum- 
nar phase exists. If r„ > T/v, then the sliding columnar 
phase only exists up to length scales of order the KT 
correlation length ^at — dexp{b/\T — T/vj-'^/^), where b is 
a nonuniversal constant. Beyond ^at, the phase will be 
a lamellar nematic. The possibility of T^ < Tn will be 
addressed elsewhere [g[ . 

In our discussion to this point, we have assumed that 
lipid bilayers define absolutely rigid 2D galleries for the 
DNA. In reality, the lipid bilayers fluctuate. To treat 
these fluctuations, we should introduce a height vari- 
able /i" (r) for membrane n and a y-displacement variable 
Uy{r) for the DNA lattice in the nth gallery. Clearly the 
variables are coupled, and at long wavelengths, we can 
assume they are locked together. Thus, we can construct 
a continuum elastic hamiltonian in terms of u(x) = Mz(x) 
and Uy{x): 



riyz — 2 



2By'u„_u„ 



-B'uly 



+KyM^v? + 2Kyjd,d.uyf + Ky^uyf]. 



The DNA lattice introduces a preferred direction in the 
xz plane that causes the bend elastic constant tensor ii'*' 
for Uy to be anisotropic. The statistical properties of u^, 
including the correlation function S'(q), predicted by this 
Hamiltonian are essentially identical to those predicted 
by the simpler rigid-layer Hamiltonian of Eq. (p. The 
Grinstein-Pelcovits renormalization of B^ , K^^ and K^y 
are identical to those of B, K, and Ky. We can use Hyz 
to calculate smectic correlations in the lipid bilayers de- 
termined by Sh{x) = (e^9o[«„(x)-«„(o)]^^ .^^j^gj.g ^^ ^ 27r/a. 

X-ray diffraction measures Sh{q), the Fourier transform 
of Shlji.), as well as S'(q). The most noticeable differ- 
ence between this function and that of standard lamellar 
phases arises from the anisotropy in the bending modu- 
lus. There are other differences arising from the u^z — Uyy 
coupling that will be discussed in more detail in an up- 
coming publication. 

The intrinsic chirality of DNA molecules leads to in- 
teractions that will cause the direction of DNA lattices 
to rotate, like the director in a cholesteric phase, about 
an axis perpendicular to the layers in both the sliding 
columnar and lamellar nematic phases. We expect, how- 
ever, these interactions to be small and resultant pitches 
to be very large because the lipid bilayer prevents close 
approach of DNA molecules. 

In this letter, we have introduced and investigated 
some of the properties of a new phase of matter that may 
exist in highly anisotropic, rotationally invariant colum- 
nar systems such as DNA-cationic-lipid complexes. It 
would be interesting to find evidence for this phase in the 
DNA complexes studied in Ref . 0] . The material in these 



experiments consisted of micron-size spherulites com- 
posed of apparently randomly oriented domains of length 
of order L = 500A on a side of DNA-lipid-complex. The 
x-ray data are consistent with nearly independent har- 
monic smectic layers with some exponentially attenuated 
positional correlation between layers. A nonvanishing 
Vu will lead to exponential correlations between layers 
even in the harmonic regime. 500A corresponds approx- 
imately to the nonlinear 2D length l^- Thus, it is not 
clear whether the domain size is set by actual distortions 
of the lipid layers or by the breakdown of harmonic be- 
havior at length scales longer than l^ or l^- It will be 
necessary to prepare monodomain samples to determine 
whether there is a crossover to anharmonic 2D behavior 
or to that of the 3-D sliding phase proposed in this paper. 
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the importance of dislocation interactions in the sliding 
phase, and to Robijn Bruinsma and Randy Kaniien for 
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the NSF under grant No. DMR94-23114. 

Note added Results essentially identical to those re- 
ported here were obtained independently by Golubovic 
and Golubovic as reported in Ref. (Tsl. 
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